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Yield of direct photons in Pb+Pb collisions at SpS and LHC energy is evaluated with emphasis on 
estimate of possible uncertainty. Possibility of experimental observation of direct photons at LHC 
is discussed. Predictions of several models at SpS energy are compared with experimental data. 



Photon emission is one of the most actively used observables both in present and developed heavy ion experiments. 
The reason is that photons, in contrast to hadronic probes, have extremely large free path length in the hot matter 
^\^ , and escape from it without rescattering. This makes them a unique probe of the initial stage of collision, and the inner 
0^ ' part of the created hot matter. In particular, direct photons are considered as one of the most promising signatures 
^""^ , of quark-gluon plasma. 

P^ ' As well known, photons, emitted in the heavy ion collision, can be divided onto direct and decay ones. The direct 

^ , photons are emitted as a result of rescatterings of charged particles in the hot matter. The decay photons are the 

' ■ photons, originated due to decays of final hadrons, mainly 7r° and 77 mesons. In addition to this, one can subdivide 

^^ , direct photons on prompt and thermal. Prompt photons are emitted in the very beginning of the collision in scatterings 

" ■ of charged constituents of target and projectile. They have approximately the same energy distribution as emitting 

I , particles - power one - and thus contribute mainly into the hard part of the spectrum. Thermal photons are emitted 

^ ' on the later stages of collision, when local thermodynamic equilibrium is almost reached. They have approximately 

l/~) exponential spectrum. Decay photons can be subdivided on photons, originated from decay of hadrons, emitted from 

CN QGP surface [|] and from hadronic gas, and therefore bring out information about two different phases of the hadron 

matter. 

Despite a large number of papers devoted to photon emission in heavy ion collisions g- |2g], it is not quite clear, 
^^ , how predictions of various models differ from each other, how strong is sensitivity of these predictions to the variations 
^^ ' of model parameters, approximations, and assumptions of the models, and how reasonable are attempts to find photon 
(-H signatures of QGP. These points are not discussed in the literature detailed enough to estimate simultaneously all 

O ' important contributions to photon spectra in wide range of photon momentum. In this paper we made an attempt 
' , to compare predictions of various models in the common basis and estimate sensitivity of the predictions to the 
bj ■ variations of the model parameters. Thus, basing on SpS data, we estimate possible uncertainties of predictions of 
(~| , direct and decay photon yield in Pb+Pb collisions at LHC. Having in mind present and developing experiments we 
L| restrict ourselves by Ph nuclei. 
. ^H : The plan of this paper is following. In the first part we consider yield of prompt photons and estimate uncertainty 

/^ • in their evaluations. In the second part we consider evaluation of emission rates of thermal photons from QGP and 
hadronic gas with and without chemical equilibrium. Having emission rates we evaluate yield of thermal photons 
accounting dynamics of the system within various hydrodynamic models and explore sensitivity to the model approx- 
imations. Then we compare predictions of cascade and hydrodynamic models at SpS energies. In the third part we 
consider decay photons and evaluate ratio Direct/Decay photons for various hydrodynamic and cascade models. Our 
results are summarized in the conclusion. 

I. PROMPT PHOTONS. 

Prompt photons are emitted in first interactions of charged constituents of target and projectile. Prompt photons 
refiect dynamics of the very beginning of collision, and can be used to study structure functions of colliding nuclei. 
To evaluate yield of prompt photons one should convolute momentum distributions of incident partons (structure 
functions of colliding nuclei) F with cross sections of elementary collisions qg — * 75, qq — > 73, bremsstrahlung etc. 

' a.h ' 

Due to the comparatively high momentum transfer, these cross-sections can be evaluated within perturbative QCD 
with reasonable accuracy. The main uncertainty in the prompt photon yield comes from structure functions: emission 



^ 



of photon even with as high energy as ii^^ ^ 5 GeV in Pb+Pb coUision at ^/s — 6300 A ■ GeV is determined by 
structure functions in the low x region [x ~ E^j 'Js ~ 0.001), where they are not known well yet. In addition, one 
should take into account modifications of structure functions in nuclei, what further increase uncertainty in the yield 
of prompt photons. 

This uncertainty is demonstrated on the fig. H, where several evaluations with different structure functions are shown. 
Dotted line corresponds to the predictions of Alam et al. Q. They used Duke and Owens, set I |Q parametrization 
of structure functions without accounting any modifications in nuclei. However, since this paper a new data on the 
low X behavior of structure functions was obtained, and in the paper ||l^, the yield of direct photons was evaluated 
with new parametrization (Gluck, Reya and Vogt pi|) of structure functions. This prediction is shown by solid line. 
As one can see, various parametrizations results in difference more then order of magnitude. In addition to this, 
structure functions in nuclei significantly differ from ones of single nucleon: depending on the value of x one finds 
shadowing, antishadowing, EMC effect, Fermi motion etc. As far as we are interesting in low x region, the shadowing 
is the most important for us. In the paper |12| Gluck-Reya-Vogt structure functions was modified to include this 
effect. This result is shown by dashed line. 




FIG. 1. Yield of prompt photons in Pb + Pb collision at LHC energy, evaluated in [n,P_2| with structure functions of Duke and 
Owens, set I - dotted line, with Gluck, Reya and Vogt parametrization of structure functions without accounting of shadowing 
(solid line) and with shadowing (dashed line). 



So that, using of different parametrizations of structure functions results in difference in the yield of prompt 
photons up to order of magnitude, while taking into account modifications of structure functions in nuclei varies 
it approximately in three times. On this background one can neglect uncertainties arising from higher order QCD 
corrections, which (estimated, e.g., from pp data) do not exceed a factor 2. 



II. THERMAL PHOTONS 



Evaluating yield of thermal photons, one assumes, that they are emitted from a medium at local thermal equilibrium. 
So that one can average individual collisions, producing photons, over time and volume and evaluate average number 
of photons, emitted from unit volume per unit time - emission rate R. In the case of perfect thermal equilibrium 
emission, the rate depends only on the temperature of the medium, however, in the absence of the chemical equilibrium 
it depends also on fugacities of quarks and gluons in QGP and hadrons in hadronic gas. When emission rates from 
QGP and hadronic gas are evaluated, one integrates these rates over space-time volume, occupied by the system, and 
finds total yield of thermal photons. 



A. Emission rate from Quarlt-Gluon Plasma 

Thermal photons are emitted from QGP due to reactions: gq — > jq (Compton scattering), qq —^ .97 (annihilation), 
qq — > ^qq (bremsstrahlung) and other processes with higher order in a = 1/137 and as ~ 0.4 - the strong coupling 
constant. To begin with, let us consider the simplest two-particle reactions: Compton scattering and annihilation. 
In this case, emission rate is described by the following formula, practically the same, as (|l|), but with thermal 
distributions instead of structure functions, Pauli blocking and Bose enhancement for final particles: 
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where sum is taken over all possible (two-particle) reactions, f{p) - thermal distribution functions, |Mp - squared 
matrix element of the corresponding reaction. 

Using \M\ evaluated at lowest order both in a and as, and taking thermal distributions, one obtains expression for 
emission rate 01: 
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This expression diverges for massless quarks m —f 0. This is reflection of the well known Coulomb divergency. 
To screen it one has to take into account medium effects. This was done in p3| using Braaten-Pisarski technique 
of hard thermal loops. As a result of this calculations current mass of quarks m was changed to screening mass 



= g^T^/2.912 {as = g^/^ir) what results in the rate: 
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Emission rate, evaluated at T = 0.2 GeV, as = 0.4 for plasma, consisting of two flavor quarks and gluons is shown 
on the fig. by dashed line. To perform these calculations analytically, Boltzmann distributions was used instead of 
Fermi and Bose ones. In case of Bose and Fermi distributions one can perforin evaluations numerically. It was shown 
in ||l^ that one can fit this result by function (||) with argument of In changed to 2.912E /g'^T + 1. Corresponding 
curve is shown on the fig. |2| by solid line, marked by rectangles. 

Let us now return to the more complicated processes - bremsstrahlung and annihilation with scattering on third 
quark. These processes arises in the next order of as, compared to Compton scattering and annihilation. It is well 
known, that in vacuum the bremsstrahlung contribute mainly into the soft part of spectrum. So one usually neglects 
its contribution into thermal photon emission. However, recently it was shown |14| , that in thermal QCD these 
processes appear to be important and even dominate over compton scattering and annihilation even in hard part of 
spectrum. In this paper in framework of hard thermal loop technique an expression was obtained for emission rate of 
hard {Ej > T) photons due to bremsstrahlung: 
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and due to annihilation with scattering on third quark: 
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where, for two quark flavors constants Jt ~ 4.45 and J^ « —4.26. It is interesting, that this contributions have the 
same order in Ug as two-particle reaction. This takes place because of the presence of collinear divergency in matrix 
element of reaction, which is screened by quark mass. This results in enhancement factor ^ i/fn^ ~ 1/ff^- Sum of 
these two contributions into thermal photon emission rate is shown on fig. y by thick solid line. As one can see, 
taking into account of these processes increases significantly (up to order of magnitude) emission rate. 



JO 



^10 

CD 

O 
^x 10 

TJ 






10 r 



10 



10 r 





^\^ : 




A 


''\^ 






" ^^^xKX 




6 


" ^^^K^ 




7 


- i i ^%;^ 


x\ 




^ : : i ^ 


xx^ 


9 


1,1,1 


T^ 



0.5 



2.5 _, ,3 

E, (GeV) 



FIG. 2. Emission rates of thermal photons from QGP at T = 200 MeV. Thick sohd hne - emission rate with accounting of 
bremsstrahlung and annihilation with rescattering on third particle. Solid line with rectangles and dashed line correspond to 
emission rates from Compton scattering and annihilation with Bose- Fermi and Boltzmann thermal distributions correspondingly. 
Dotted line and thin solid line - chemically nonequilibrated QGP with Xq = 0.5, Ag = 1 and Xq — 0.5, Xg — 0.5 correspondingly, 
only Compton scattering and annihilation are taken into account. 



If QGP is formed in the nucleus-nucleus collision, then the number of quarks and gluons in it may be much lower, 
than in the plasma at perfect chemical equilibrium. This can strongly affect on the emission rate Pq-p?7| . The simplest 
way to account undersaturation of quarks in the QGP is just to put coefficient A, (fugacity of quarks) as a coefficient 
in the emission rate. More advanced possibility is in addition to this to evaluate screening masses, depending on 
quark and gluon fugacities. However, in the case of absence of chemical equilibrium one can not be sure that there 
will be perfect screening of the column divergency as in the equilibrated case. In the paper |19| it was shown, that 
even in the case of absence of chemical equilibrium divergencies cancels each other in the emission rate. In particular 
case of distributions of quarks and gluons: 
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it was shown, that screening mass is 
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and emission rate is given by formula 
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of Rieman C-function (C'(2) = -2.404). 

In the case Ag = 1 and Xg = 1 this formula coincides with (0). Emission rates, evaluated with this formula are 
shown on fig. by dotted line (A^ = 0.5, \g = 1) and solid line (A, — 0.5, \g — 0.5). As one can see, emission rate is 
very sensitive to the undersaturation of charged constituents, while undersaturations of gluons slightly vary emission 
rates. 

There are three main sources of uncertainties of emission rate. First - contributions from higher order in as 
processes. As we have seen, taking into account bremsstrahlung and annihilation with scattering on third particle in 
addition to Compton scattering and annihilation result in increasing of emission rate up to order of magnitude. So 
one can not be sure that inclusion of processes of the next order in as does not result in similar changes of emission 
rate. Second - all considered evaluations are performed in the framework of Braatcn - Pisarski Hard Thermal Loop 
technique, applicable for the case 5 <C 1. To evaluate yield of thermal photons in the heavy ion collisions we have 
to extend these results for the case as ~ 0.4 e.g. g ~ 2. Third - in our evaluations we consider region Ej 3> T - 
the right tail of the thermal distribution. It is well known, that thermal distribution on the tails forms much more 
slower, than in the region _E ~ T, where 1-2 collisions are sufficient for establishing of thermal equilibrium. So, if 
temperature is rapidly changes with time, then long tail is not able to form, and it is the region, where this effect is 
important. In contrast to the first two sources, the last effect results in changing of the shape of the distribution and 
can even result in nonexponential distribution. 



B. Emission rate from hadronic gas 



Thermal photons can be emitted from hadronic gas as a result of reactions tttt — > p7, np — > ttj, p — > 7r7r7, tttt — > 7/7 
etc. In the paper |13] it was shown, that processes Tip — > 7r7 plays the dominant role in emission of photons with 
energy larger than ^ 0.7 GeV . In addition, it was shown, that photon emission of the hadronic matter is very similar 
to the photon emission of QGP (including only Compton scattering and annihilation) at the same temperature. So 
one can parametrize photon emission rate from hadronic gas in the form (HV Emission rate of hadronic gas, evaluated 
at r = 180 MeV in accordance with this parametrization is shown on figT^ by solid line. 
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FIG. 3. Emission rates of thermal photons from hadronic gas at T = 180 MeV , evaluated with emission rate i}^ (solid line), 
with accounting of ai resonance (H) (dotted line), and with accounting of in- medium effects (dashed line). 



Later it was shown [ |l8| , that contribution of Ai resonance into process np ^r Ai ^> ttj is important: it increase 
emission rate approximately twice. There was proposed following parametrization of emission rate with accounting 



of contribution of Ai resonance: 

S-^^= 2.4 T^i^exp (-1/(1.35 T^)0-^^-^/T) (fm-^GeV-^) (9) 

Considered emission rates are evaluated with vacuum masses and widths of hadrons. However, these properties 
change in medium, so that evaluated emission rates should change too. In the paper |19{ an attempt was made to 
evaluate photon emission rate with accounting of such effects. Effective Lagrangian approach was used to evaluate 
in-medium masses and widths. It was shown, that due to decreasing of vector meson masses yield of thermal photons 
increases. These predictions for hadronic gas at T = 180 MeV are shown on the fig. by dashed line. 

From one hand, situations with emission rate from hadronic gas is almost so shaky as with emission rate from QGP 
- as far as there is no theory of strong couplings. From the other hand, one can normalize coupling constants and 
other hadronic properties to experiment as it was done in [|l8| , what can not be done in the case of QGP. In addition 
to this, evaluated emission rates suppose chemical equilibrium. However, it may not be established due to various 
reasons: fashion of hadronization of QGP, fast cooling of hadronic gas, fly out of hadrons from QGP surface ||l[ and 
hadronic gas, etc. As it was shown e.g. in [|l3| , main contribution to emission rate comes from pion-pion and pion-rho 
scatterings. So, the most important for us is oversaturation of pions and rho due to fast cooling. One can estimate 
this oversaturation as ratio N{Tc)/N{Tf) ^ (Tc/Tf)^, where Tc and Tf - transition and freeze-out temperatures 
correspondingly. So one could expect, that emission rate increase in the same ratio. So we estimate uncertainty in 
emission rate from hadronic gas within factor ~ 3 due to in-medium effects and factor ~ 3 due to oversaturation of 
pions and rho in hadronic gas. 

C. Yields of thermal photons 

To evaluate yield of thermal photons in the nucleus-nucleus collision one should convolute emission rate with 
space-time evolution of the system: 
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Below we consider several models of evolution of hot matter and evaluate yield of thermal photons in these models. 

First of all, on the example of very simple 1-1-1 (1 space-|- time) Bjorken hydrodynamic model pfl] wc show main 
phenomena: dependence on initial conditions, transition temperature etc., and then compare predictions of Bjorken 
model with 1-f-l Landau hydrodynamics and more advanced model such as (2-1-1) Bjorken model. 

In the 1+1 Bjorken hydrodynamic model of evolution of heavy ion collision it is assumed, that hot matter, created 
in the beginning of collision, expands only longitudinally (along collision axis) and longitudinal velocity depends on 
longitudinal coordinate in the following way: 
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This parametrization of velocity result in flat rapidity distribution of fluid elements, and as a result - of final 
hadrons. So, this model is well applicable only in the midrapidity region. Below we restrict ourselves by midrapidity 
region z — 0, t — t, because this region is most interesting from point of view of observation of hot hadronic matter. 

As a result of such expansion, the volume of the hot matter increase with time in the following way: 

F = 2 7ri?2 t, 

where R - radius of colliding nuclei. If one assumes now, that entropy is conserved during all the expansion, one can 
relate temperature and time: 

or 
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where T - temperature, t - time, g - degeneracy, dN/dy - multiplicity of final (massless) hadrons (bosons) and 
C(3) = 1.202. In the simplest version of Bjorken hydrodynamics we have ggcp = 16 + |24 for degeneracy of QGP 
and giiadr — 3 for hadronic gas with dominating pions. 

Now, choosing from some considerations e.g. from cascade model multiplicity dN/dy and initial time (initial 
temperature) - the time (temperature), at which we already can assume thermal equilibrium, we find initial conditions. 
Then, fixing transition temperature T^ and freeze-out temperature Tf, we completely define this model. 

Having emission rate of photons (e.g. (0)) and description of evolution, one can evaluate the total number of 
thermal photons as function of model parameters: 
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As one can easily see, the total number of thermal photons slightly depends on the initial temperature, but strongly 
- on freeze-out one. This is due to emission of very soft photons during prolongated evolution in the case of lower Tf. 
However, we are interested in the hard part of the spectra, so that this quantity is not appropriate for our purposes. 
More informative would be pt distribution of thermal photons. 

Let us now fix basic set of parameters of the model: 



dN/dy = 12000, T,n = 1 GeV, T^ = 160 MeV, Tf = 100 Mel^, 
and look, how pt distribution of direct photons depends on variation of model parameters. 
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FIG. 4. Yield of thermal photons in Ph -\- Pb collision at LHC energy within l-fl Bjorken hydrodynamics. Left plot: yield 
evaluated using emission rate (^j both for QGP and hadronic phase and parameters ( |l2|) for dN/dy — 12000 (solid line), 
dN/dy = 9000 (dotted line) and dN/dy — 6000 (dashed line). Right plot: Yield, evaluated with parameters (hd) and initial 
temperatures: Tin ~ 1 GeV (solid lines), Tin ~ 500 AleV (dashed lines) and Ti„ — 300 MeV (dotted lines). Upper three curves 
evaluated using emission rates (m),(p|) for QGP phase, and (Bl) for hadronic phase, bottom three curves - using emission rate 
(W) both for QGP and hadronic phase. 



To begin with, let us vary multiplicity - see fig. 0. With fixed initial temperature variation of dN/dy means variation 
of initial time, and thus all times of the evolution: all times change in the same scale. As one can expect from ([ll]), 
the number of photons varies as squared multiplicity, while the shape of distribution does not changes. However, this 
effect holds only for l-fl Bjorken hydrodynamics. In more realistic models the dependence of dN'' /dy on the total 
multiplicity is weaker. For example in the 2-1-1 bjorken hydrodynamics N'' ~ (dN/dyY'^ pl[. This can be understood 



in the following way. The number of thermal photons is proportional to the number of charged particles times the 
number of collisions, suffered by each particle: 

N^ - Nchiy. 

For the case of long lived system each particle can collide to all other particles, so i^ ~ Nch- However, if system 
expands, or has final dimensions comparable or even less then free path length, then the number of collisions 

D 

where D - dimensions of the system. In the case of 1+1 Bjorkcn hydrodynamics increasing of dN/dy result in 
increasing of longitudinal dimensions of the system at fixed transverse sizes, so I? ^ dN/dy 

N^ - (dN/dyf. 

If we allow transverse expansion, then D ~ {dN/dyY^^ and 

N^ ^ {dN/dyf'\ 

result very close to obtained in 2+1 Bjorken hydrodynamics. 

As a next step we vary initial temperature - see fig. 0, right plot. With fixed dN'" /dy this means that we cut 
off initial stages of evolution: in the case Tin — hQOMeV we skip time r — 0.15 — > Ifm/c with respect to the 
case Tin = IGeV and in the case Tin = 300 MeV^ - r = 0.15 — > 4.6 /?ti/c. So, from fig. |j one can find what part 
of spectrum is formed at each stage of collision. We obtain very natural result: the harder part of spectrum is 
populated due to emission from the hottest stage of collision. In addition, on this figure we compare yields of thermal 
photons, evaluated using simplest emission rate (0) (three lowest curves) and emission rates, accounting second order 
QCD corrections (H,H) and ai contribution (m (three upper curves). In the latter case yield of thermal photons is 
approximately order of magnitude larger. 
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FIG. 5. Left plot: The same as fig. 0, but with different transition temperatures Tc — 200 MeV (dashed line), Tc = 180 MeV 
(dotted line) and Tc — 160 MeV (solid line). Right plot: thermal and decay photons evaluated within 1+1 Bjorken hydrody- 
namics with different freeze-out temperatures: T/ = 100 MeV (solid line), T freeze = 120 MeV (dashed line) and Tc = 140 MeV 
(dotted line). In all cases thermal photons are evaluated with emission rate 



Next step is variation of the transition temperature Tc - see fig. |5| left plot. As we see below, mixed phase lives 
much more longer, than QGP phase, and thus a lot of thermal photons are emitted with T = Tc, increasing yield of 
medium energy direct photons with increasing of Tc- Lifetime of mixed phase is determined by ratio of degeneracies 
of QGP phase and hadronic phase. The simplest Equation of State (EOS) of hadronic gas corresponds to ideal gas of 
massless pions, so that this ratio is 37/3 ~ 12 and mixed phase lives a very long time. If one considers hadronic gas, 
consisting from all known hadrons with masses below 2GeV, then this ratio becomes 37/12 ^ 3 and time of life of 
mixed phase is significantly reduced, reducing yield of medium energy photons. 

And finally we vary freeze-out temperature Tj. Its variation affects both pt distributions of thermal photons and 
Pt distributions of final hadrons, and thus decay photons - so that we show both decay and direct photons on the 
fig. ra, right plot. From one hand, decreasing of Tf result in increasing of time of life of hadronic gas and thus yield 
of soft thermal photons, from the other hand it result in decreasing of the temperature of final hadrons, so that ratio 
thermal / decay increase. 

Comparing figs. |^^ one can conclude, that within 1+1 Bjorken hydrodynamic model hard part of spectrum 
of thermal photons pt > 2GeV/c is governed by emission from the hottest region (QGP), variation of transition 
temperature vary pt distribution at 1 < pt < 2.5GeV/c, while the softest part of the spectrum pt < 0.5 GeV/c 
depends on the freeze-out temperature. 

It is remarkable (see fig. |5|, right plot), that within 1+1 Bjorken hydrodynamics direct photons dominate over decay 
ones practically at all energies. We obtain this unphysical result because we neglect transverse expansion. In the case 
of thermal photons absence of transverse expansion result in extremely large times of evolution: e.g. for our set of 



parameters (12) we find 

To = 0.15 fm/c, Tqgp = 12 fm/c, r,„„ed = 137/to/c, Tf = 560 fm/c, (13) 

where tq - time of thernialization, Tqgp - time of disappearing of QGP, Tmixed ~ time of disappearing of mixed phase, 
Tf - time of freeze-out. So that, emission of photons from the mixed phase and especially from the hadronic phase is 
strongly overestimated in comparison with the models with transverse expansion. 

In the case of decay photons, absence of transverse expansion result in lower effective temperature, and smaller 
contribution with respect to thermal photons. 

However, this docs not mean, that 1+1 Bjorken model is useless: as we see below, it makes reasonable predictions 
for hard part pt > 2 GeV/ c of the spectrum of thermal photons, emitted on the stage, where transverse expansion is 
not very important. 

Bjorken model of evolution is essentially nonhydrodynamic, i.e. it is not solution of the hydrodynamic equations. 
So, it is interesting to compare its predictions with really hydrodynamic model e.g. 1+1 Landau hydrodynamics 
p2| . Within this model it is assumed, that just after penetration of colliding nuclei through each other a thermalized 
matter is created, having form of Lorentz compressed disk. The part of the total energy, deposed in this disk is 
model parameter - coefficient of inelasticity K . Being formed, disk expands along beam axis in accordance with one 
dimensional hydrodynamic equations. Initially disk has size along beam axis: 

4m^^ (14) 

where nip is the proton mass, K ~ coefficient of inelasticity K ^ 0.5, R - radius of nucleus and cj. - speed of 
sound in the created matter. For Ph + Pb collision {R = 6.5 /m) at LHC energy (-^s = 6300 A • GeV) one finds 
(/^ = 0.5,c2 = l/3): 

Aw 0.0025 /m. 

Initial temperature can be evaluated from the condition 

27ri?2A • £ = Ky/1 ■ A, 

where £-energy density. Assuming QGP consisting of two quark flavors and gluons, one finds for Ph + Ph collisions 
at LHC energy: 

T,^ = 5.lGeV. 

Result of evaluation of the yield of thermal photons from QGP phase only {Tc = 220 MeV) within Landau model 
taken from p3] is presented on the fig. H by solid line. Dashed line corresponds to the evaluation within 1+1 Bjorken 



model with the same initial conditions: Tin = 5.1 GeV, T^ = 220 MeV and Ti„ = 0.0025 /m/c. This conditions 
correspond to multiplicity dN^ /dy = 5300. As one can see, lines goes rather close to each other, though there is some 
difference due to different rapidity distributions and, as a consequence, different dependence T{t). Scaling expansion 
result in faster cooling in the beginning and slower - on the later stages. Nevertheless, this result demonstrate, that 
use of scaling longitudinal expansion instead of hydrodynamic one is not critical for our purposes. 
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FIG. 6. Yield of thermal photons, emitted from QGP phase in Pb + Pb collision at LHC energy, evaluated within Landau 
hydrodynamic model (solid line), and 1+1 Bjorken hydrodynamics with To = 5.1 GeV, Tc = 220 MeV and Ti„ = 0.0025 /m/c 
(dashed line). 

Taking into account transverse expansion alters significantly all evolution of the system, and thus spectra of the 
thermal photons. To include transverse expansion we consider 2+1 Bjorken hydrodynamics. This means, that we 
assume Bjorken scaling expansion along beam axis, while for expansion in transverse direction we solve hydrodynamic 
equations. So, we parametrize 4-velocity in the form: 



-r -r ^ ' 



2\-2 



and solve hydrodynamic equations 



a^T^-^^O, T^''{x)^{e{x)+p{x))-u^^ 



i'^ — ^/^^ 



P{x), 



(15) 



where T^'^ - energy-momentum tensor of ideal fluid, e{x) and p{x) - local energy density and pressure correspondingly, 
yM _ 4-velocity and g^'^ metric tensor {g = diag{l, — 1,— 1,— 1)). In the case of cylindrically symmetrical transverse 
expansion and longitudinal scaling expansion one can rewrite dl3) in the form 

dr (t^w) + dr (7^t;w) + j'^wt^^ + ■f'^vwr^^ — drP = 

dr [j'^vw) + dr (7^w^w) + j'^vwT^^ + j'^v'^wr^^ + drP = 0, 

where w = e + p - enthalpy density. We solve this equations numerically. 
As in the one dimensional case we fix basic set of parameters 



(16) 



dN/dy = 12000, T,„ = 1 GeV, T^ = 160 MeF, Tf = UQMeV. 



(17) 



In contrast to the 1+1 dimensional case, where all time scales are given by initial time (at fixed initial temperature), 
in the 2+1 dimensional case we have two completely different time scales: one of longitudinal expansion - the same 
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as in 1+1 dimensional case and one of transverse expansion, determined by velocity of sound and dimensions of 
the system. So that, picture of expansion become much more complicated. Typical evolution of hot matter in 2+1 
dimensional case is shown on the fig. - this is result of evaluations with parameters, listed in ( p7| ) (left plot) and 
dN/dy = 6000, Ti„ = 1 GeV, Tc = 160 MeV, Tf = 140 MeV (right plot). Initially system has zero transverse collective 
velocity, than due to acceleration in the QGP phase matter possesses velocity more then 0.7 • c. In the mixed phase 
pressure gradient is absent and velocities are constant (velocity levels are straight), and after the end of the mixed 
phase during hadronic gas expansion matter begins to accelerate further (slopes of velocity levels changed). 
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FIG. 7. Evolution of the hot matter, calculated at LHC energy for dNjdy = 12000 left plot, and dN/dy = 6000 right plot. 
Constant velocity levels (thin lines) are drown at 0.1,0.2,... and constant energy density levels, are drawn at the boundary of 
QGP phase (thick solid line), at the boundary of mixed phase (thick dashed line) and freeze-out surface (thick dotted line). 



In contrast to one dimensional case, times of life of mixed and hadronic phases is much smaller: 



To = 0.11 fm/c, Tqgp = 10/to/c, Trmxe, 

in the case of dN/dy — 12000, and 



d = 25 fm/c, Tf = 30 fm/c. 



To = 0.06 fm,/c, 



i fm/c, 



ed = 22 fm/c, Tf = 25 fm/c. 



(18) 



(19) 



in the case of dN/dy = 6000. 

In the case of larger multiplicity initial time is twice larger, so that in the beginning stage system cools twice slower, 
and more energy is spent to the developing of the radial velocity. Thus we find higher transverse collective velocity in 
the case of higher multiplicity. However, later stages (i > 5 /m/c) are governed mainly by radial expansion, so that 
life times of QGP, mixed phase and hadronic gas are very similar on both cases. 

There are two important consequences of introduction of transverse expansion. First, - times of life are changed: 
time of life of QGP is changed slightly, but time of life of mixed phase and hadronic gas are reduced more then order 
of magnitude compared to 1+1 dimensional case. Second, ~ hot matter poses significant radial velocity, which affect 
spectrum of thermal photons and even in the more degree ~ spectrum of decay photons. 

As a result of producing of the large radial velocity the thermal photons, emitted from hadronic gas contribute 
into the hard part of the spectrum of direct photons - see fig. H. On these figures contributions of different times 
into the total yield of thermal photons are shown. Similar to one dimensional case, photons emitted at first several 
fm/c populate hard part of spectrum. However, due to high collective velocities contributions of later (cooler) stages 
into hard region are comparable with emission from the hottest stage. This effect is better pronounced in the case of 
higher multiplicity {dN/dy = 12000), where collective velocity can reach Vr ~ 0.9 and effective temperature 
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T, 



eff 



Tf 



0.14 Gey -4.4- 0.6 GeF 



(20) 



is comparable even with initial temperature. 
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FIG. 8. Contributions of different times (from down to up 0-1, 0-5, 0-10, 0-15, 0-20, 0-25 fm/c, total) into the total thermal 
photon yield, evaluated at dN/dy = 12000 (left) and dN/dy = 6000 (right). 



Similar to 1-1-1 dimensional case we investigate dependence on the parameters of the model. In particular on initial, 
transition and freeze-out temperatures. Yields of thermal and decay photons, evaluated with three different initial 
temperatures {Tin = IGeV, 500 MeV and 300 Afe^), and fixed multiplicity, transition and freeze-out temperatures 
(O) are shown on the fig. 0, left plot. If one compare spectrum of thermal photons, evaluated within 1-1-1 (fig. 0, 
right plot) and 2+1 (fig. 0) Bjorken hydrodynamics, one finds, that in the case of Tin = 1 GeV hard parts of spectra 
of thermal photons coincide. So one can conclude, that both in 1+2 and 1+1 cases the system on initial stages 
cools predominantly due to longitudinal expansion. For lower initial temperatures Ti„ — 500 AleV and especially 
Tin = 300 Mey thermal photon spectra evaluated within 2+1 hydrodynamics are harder. This is evidence, that 
emission from later stages with lower temperature, but higher radial velocity and thus high effective temperature (pfl) 
dominate over emission from QGP. 

Sensitivity of the yields to the variations of transition and freeze-out temperatures are shown on the fig. 0, right 
plot. Evaluations with basic set of parameters (O) (solid lines), higher transition temperature Tc — 200 MeV (dotted 
lines), and lower freeze-out temperature Tf = 120 MeV (dashed lines) are shown. In all three cases decay photons 
goes very close to each other. The same are thermal photons, however decreasing of Tf slightly increase yield of 
intermediate thermal photons. 

In all previous calculations we used the simplest equation of state (EoS) - EoS with first order phase transition 
with ideal gas of quarks and giuons (degeneracy qqgp — 37) in one phase and ideal gas of massless pions (degeneracy 
ghadr = 3) in thc othcr phase. Such EoS includes large jump in of specific entropy, and thus result in long living 
mixed phase. Besides this, such EoS leads to maximal possible velocity of sound v^ — •\/l/3 both in QGP phase and 
hadronic phase and thus maximal possible radial velocity at fixed value and spatial distribution of initial energy. If 
one considers massive hadrons instead of massless, and interacting quarks and giuons instead of free ones, then one 
finds Vc < y^l/3. In addition, this inclusion of heavier hadrons and resonances in hadronic gas significantly increase 
its degeneracy and reduce jump of entropy. These effects result in decreasing of the yield of thermal photons in 
intermediate region. 
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FIG. 9. Left plot: yield of thermal and decay photons, evaluated within 2+1 hydrodynamics with set the of parameters ( |l7| ) 
and initial temperature Tin = IGeV (solid lines), Tin = 500 MeV (dotted lines) and Tn = 300 Af eV^ (dashed lines). Right 
plot: predictions with the basic set of parameters nVJn (solid lines), variation of transition temperature Tc = 200 MeV (dotted 
lines) and freeze-out temperature Tf = 120 Af eV^ (dashed lines). In all evaluations emission rates (0) (n) and (bf) are used. 

Another phenomenon perhaps important in evaluation of direct photon emission in ultrarelativistic heavy ion col- 
lision is absence of chemical equilibrium in the QGP phase. It was argued |15| , p3| -]2^] that in the very beginning of 
collision quarks and gluons are strongly undersaturated (approximately 0.1 of equilibrium multiplicity), and subse- 
quently their number is increased due to reactions gg <-^ ggg, gg <->■ gg etc. Because of comparatively small coupling 
constant, chemical equilibration is reached rather slowly - see fig. lOl left plot, with results from ref. pa. In these 
evaluations initial conditions are taken from HIJING event generator, and then Boltzmann equation was solved with 
collisional term including all possible QCD binary interactions as well as reaction gg <-> ggg, included to describe 
establishing of chemical equilibrium. Solid lines on the fig. nu correspond to the evaluations with fixed Og — 0.3, and 
dashed lines - to evaluations with running coupling constant 



as{<E>) 



47r 



9\n(^<E>^ /Al^^ 



where <E> - average energy of partons evolved. As one can see from this figure, due to higher coupling, gluons come 
to the equilibrium faster than quarks and rate of equilibration practically independent on a^. In contrast to them 
quarks equilibrate much slower, and this process strongly depends on as- However, one should take these evaluations 
with care: cross sections, entering collisional term has Coulomb divergency. To screen them one has to introduce cuts 
on soft quark and gluon exchange from some considerations, and result should be very sensitive to the chosen values 
of these cuts. 

Yield of thermal photons emitted from QGP, evaluated within this approach with accounting of Compton scattering 
and annihilation and with time dependent a^, is shown on fig. nOl right plot, by solid line. We compare this yield with 
one, evaluated within 2-1-1 Bjorken hydrodynamics with different initial temperatures and using emission rate (Q) i.e. 
also including only Compton scattering and annihilation. We find that in the case of absence of chemical equilibrium 
yield of thermal photons decrease by lO^'^ — 10^^. 
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p, (GeV/c) 

FIG. 10. Left plot: fugacity of quarks and gluons at LHC, solid lines - fixed as, dotted lines - running a^. Right plot: thermal 
photon yield from QGP without chemical equihbrium (solid line) and from 2+1 Bjorken hydrodynamics with Ti„ = IGeV 
(dashed line) Tin = 500 MeV (dotted line) and Ti„ — 300 MeV (dash-dotted fine). Only Compton scattering and annihilation 
(emission rate (W)) are taken into account. 

However absence of chemical equilibrium does not influence significantly on the yield of thermal photons. As we 
have seen on fig. pi due to significant radial velocity contribution of the hadronic gas, it is approximately the same 
as contribution of the hot phase. So even disappearing of the emission from the hottest stage does not reduce total 
yield of thermal photons more than several times. 

D. Thermal photons from QGP and hadronic gas 

From the point of view of registration of QGP it is important to know, whether or not thermal photons from QGP 
dominate over the thermal photons from hadronic gas. From one hand, thermal photons emitted from QGP have 
highest temperature, from the other hand thermal photons from hadronic gas have higher radial collective velocity, so 
that the answer is not obvious. In the paper m it was shown, that within 2+1 Bjorken hydrodynamics and emission 
rate ( fm th e contribution from QGP is slightly smaller than contribution from hadronic gas. However, in the recent 
paper |26[| , emission of thermal photons was evaluated with a rate accounting bremsstrahlung (|5|) and annihilation 
with scattering on third particle (|g) . It was shown, that in this case even for initial temperatures Tin ~ 500 MeV 
emission from QGP dominate over emission from hadronic gas. It is interesting, that in this case energy carried out 
by thermal photons becomes comparable with energy carried by final hadrons, and evolution of the system should by 
corrected to take into account energy loss due to thermal photon emission. 

One can estimate yield of thermal photons from QGP and hadronic gas from fig. ^. Third curve from the bottom 
corresponds to the time (r < lOfm) of life of QGP phase. This curve should be compared with total photon yield. 
In the case of initial temperature Tin = 1 GeV and dN/dy — 12000 (left plot) contribution of QGP is approximately 
twice smaller than total photon yield, what means, that contributions of QGP and hadronic gas are comparable even 
a.t pt ^ 5 GeV/c. In the case dN/dy = 6000 (right plot) relative contribution of hadronic gas is smaller and QGP 
photons could dominate. 
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E. Direct photons at SpS 

Below we compare spectra of thermal photons in Pb+Pb collisions at 158 A • GeV (SpS), predicted by various 
models, with preliminary experimental results, obtained by WA98 collaboration. Our goal is not to choose model, 
which better describes experimental results or even to conclude, whether or not QGP was formed in these collisions. 
Wc just would like to demonstrate, that calculations within 2+1 Bjorken hydrodynamics can reproduce measured 
spectrum of thermal photons, providing we guessed appropriate initial conditions. As well we show, that at least at 
this energy different models give very similar predictions, so it is possible, that at LHC energy these predictions will 
be close too. 

Let us begin from 2+1 Bjorken hydrodynamics. In this case we fix Tc — 160 MeV, Tf = 140 MeV and choose initial 
conditions Tm ~ 320 MeV, Tin = 0.4 fm/c so that reproduce multiplicity dN/dy = 620 and spectrum [g^ of final 
photons. Quality of the data description can be seen on the fig. O, left plot. Having parameters of the model fixed, 
we evaluated yield of thermal photons using two emission rates of thermal photons: in the first case we use emission 
rate M) both in QGP and hadronic gas phase and in the second case - emission rate (O), (ph for QGP phase and (H) 
for hadronic gas phase. Evaluated yields of thermal photons are shown on the fig. nM, right plot, by thin dashed line 
and thin solid line for the two cases respectively. Experimental data are shown by rectangles. 

Second hydrodynamic model we consider is fully hydrodynamic (3+1 dimensional) model |^. Yield of thermal 
photons, evaluated within this model is shown by thick solid line. In these evaluations emission rate (^ for QGP 
and (Q) for hadronic gas was used. However, the difference between predictions of these models arises not from type 
of longitudinal expansion - scaling or hydrodynamic, but from initial conditions: in the 3+1 hydrodynamics they 
are evaluated in spirit of Landau hydrodynamics (shock waves propagate through the colliding nuclei and heat up 
the system) . Such description of initial stage of collision results in higher initial temperatures compared to usually 
considered in Bjorken scenarios, so one finds higher effective temperatures of thermal photons. 
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FIG. 11. Left plot: total photon yield in Pb+Pb collision at 158 A- GeV, measured by WA98 collaboration, preliminary data, 
(crosses) and evaluated within 2+1 bjorken hydrodynamics (dashed line). Right plot - yield of direct photons in this collision. 
Rectangles - preliminary WA98 results, dashed histogram - predictions of Parton Cascade Model, triangles - predictions of 
UrQMD, thick solid line - 3+1 hydrodynamics, thin solid and dashed lines - 2+1 Bjorken hydrodynamics with emission rate 
accounting of bremsstrahlung and contribution of ai resonance and without accounting of these effects correspondingly. 

To compare predictions of equilibrium (hydrodynamic) models with nonequilibrium (cascade) models we consider 
as well predictions of Partonic Cascade Model (PCM) |g8[ and Ultrarelativistic Quantum Molecular Dynamic Model 
(UrQMD) model ||2^. PCM describes nucleus- nucleus collisions in terms of parton picture of hadronic interactions. It 
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performs evolution using quantum kinetics with interactions of quarks and gluons, evaluated within pQCD. Developed 
partonic picture clasterizes to form final hadrons. Predictions of this model for Pb + Pb collisions at 160 A ■ GeV, 
taken from [g8[, are shown on fig. |ll|, right plot, by dotted histogram. In contrast to PCM, the second cascade 
model, UrQMD, deals with propagation and interactions of hadrons in hot matter. Yield of direct photons in Pb + Pb 
collisions at 160 A • GeV evaluated within this model, taken from |29], is shown on fig. O, right plot, by triangles. 

Comparing predictions of all considered models one finds that they go rather close to each other. Moreover, com- 
paring PCM and UrQMD one finds a striking coincidence of their predictions. It looks like the yield of direct photons 
is independent on the details of the interactions of particles in the hot matter. As for 2+1 Bjorken hydrodynamics, 
one can see, that yield, evaluated with accounting of bremsstrahlung and annihilation with rescattering on third 
particle gives the best description of experimental data among considered models. So we conclude, that 2+1 Bjorken 
model provides reasonable description of evolution of hot matter, such as we obtain good reproduction of the yield of 
thermal photons. 

Another important observation, we can made from fig. |ll|, is that in the range pt < 5 GeV/c predictions of all 
considered models differ less than order of magnitude. 



III. DECAY PHOTONS 



In the previous sections we evaluated yield of thermal photons. However from the experimental point of view not 
only its absolute value but also ratio of spectra of thermal and total photons is important. It is this ratio, which 
controls possibility of extraction of thermal photons from the total photon yield. To estimate total photon yield we 
have to evaluate spectrum of decay photons, i.e. to evaluate pt distributions of tt*^, tj, etc. hadrons, having channels 
of decay onto photons, and decay them into photons. To do this we use several models: 2+1 Bjorken hydrodynamic 
model, HIJING 1.36 and VENUS 4.12 event generators. 

On the fig. K2 we present results of hydrodynamic calculations of yields of decay photons for three values of initial 
temperature Tin = IGeV, Tin = 500 MeU and Tin — 300 MeV, and compare them with predictions of HIJING 1.36 
and VENUS 4.12 models. 
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FIG. 12. Total photon yield in Pb + Pb collision at LHC energy, evaluated within 2+1 Bjorken hydrodynamics with initial 
temperatures Tin = IGeV (solid line), Tn = 500 MeV (dashed line) and Tn = 300 Me V (dotted line). Crosses - predictions 
of HIJING 1.36, circles - predictions of VENUS 4.12. 

HIJING is model of heavy ion collisions, which consider initial stage of collision (jet and minijet production) on the 
basis of perturbative QCD, and describe fragmentation of created partons using Lund string model |3^ . HIJING does 
not include interactions of the secondary particles. This results in nonexponential pt distribution with excess of high 
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energy and lack of intermediate energy particles in comparison with thermal distribution. In contrast to HIJING, 
VENUS takes into account rescattering of final particles 1 33 1 and thus results in exponential spectrum, closer to the 
predictions of hydrodynamic model. However, VENUS does no include radial flow, so that spectrum predicted by 
this model is softer than one, predicted by hydrodynamics. In evaluations of the ratio thermal to decay photons we 
use predictions of spectrum of decay photons from hydrodynamic model and HIJING. 

If one uses pt distribution of decay photons, given by HIJING 1.36, (see fig. O, left plot) then one finds, that thermal 
photons dominate over decay photons in the intermediate region for all considered initial temperatures. This takes 
place because of rather special shape of pt distribution of final hadrons, predicted by HIJING. In contrast to this, if 
one evaluates both thermal and decay photons within hydrodynamic model (see fig. O, right plot), then for initial 
temperatures ^ 300 MeV there will be region, where ratio reached its maximal value and position of the maximum is 
shifted to hard part with increasing of initial temperature, while at Tin '^ 500 MeV becomes approximately constant. 
For higher initial temperatures contribution of thermal photons increase with increasing of pt ■ 

If one consider only leading order Ug reactions in QGP - Compton scattering and annihilation, then the ratio 
Thermal/Total reaches 20% in maximum even for lower considered initial temperature. If one includes higher order 
as corrections - bremsstrahlung ann annihilation with rescattering on third particle, then the ratio Thermal/Total, 
evaluated for T,„ = 300 MeV reaches 30-40%. 
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FIG. 13. Ratio yield of thermal to yield of decay photons in Pb + Pb collision at LHC energy. Left plot: Thermal photons, 
evaluated within 2+1 Bjorken hydrodynamics, are divided by decay photons, taken from HIJING. Rectangles correspond to 
thermal photons, evaluated with initial temperatures Tin — 1 GeV, diamonds - Ti„ — 500 MeV and triangles - Ti„ = 300 MeV, 
in all cases emission rate M) was used. Right plot: both thermal and decay photons are evaluated within 2+1 Bjorken hydro- 
dynamics. Solid lines correspond to initial temperature Tn — 1 GeV, dashed - Ti„ = 500 MeV and dotted - Ti„ = 300 MeV, 
lower lines evaluated with emission rate (M) , upper - with emission rates (H) , (o) and (bh . 

IV. CONCLUSION 

We have considered emission of prompt, thermal and decay photons in Pb + Pb collisions at SpS and LHC energies 
and estimated uncertainty in predictions of their yields. We find, that main uncertainty in the yield of prompt photons 
comes from the lack of precise knowledge of structure functions at small x and their modifications in nuclei. This 
leads to variations of the yield of prompt photons within factor ~ 3. 

For thermal photons situation is more complicated. The main source of uncertainty comes from impossibility to 
describe uniquely evolution of heavy ion collisions. For description of evolution the both kinetic and hydrodynamic 
models are used. We compared predictions of kinetic and hydrodynamic models at SpS energy, however for LHC 
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energy we concentrated on hydrodynamic models. We find, that hard part of the spectrum pt > 4 — 5 GeV sUghtly 
depends on the presence of the radial flow (compare figs. ^ for 1+1 expansion and ^ for 2+1). This takes place because 
this region oipt is populated due to emission from the initial stage (see fig. reftimes) when radial flow is not developed 
yet. Even if extremely high initial temperatures are created in the beginning of the collision, they are slightly affect 
on the radial flow, because as a result of longitudinal expansion the high temperatures lives very short time and can 
not accelerate considerably the matter in radial direction. This is illustrated on the fig. 0: difference of the effective 
temperatures of decay photons, evaluated at Tin — I GeV and Tin — 500 Afe^ is smaller, than difference between 
Tin = 500 MeV and Tin = 300 MeV. So, we find that for hard part of spectrum 1+1 and 2+1 hydrodynamics models 
give similar predictions. However, 2+1 Bjorken hydrodynamics is not quite two-dimensional - longitudinal expansion 
in it is scaling. Comparing predictions of 1+1 Bjorken and 1+1 Landau hydrodynamics wc show, that this is not 
important for the case of thermal photons. 

In addition to ambiguity in description of evolution, there is significant uncertainty in the emission rate: inclusion 
of the higher order a^ corrections result in increasing of it up to order of magnitude. 

We show, that hard part of thermal photon distribution is sensitive to the scenario of collision. If fast thermalization 
of colliding nuclei with high initial temperature takes place, than one finds high yield of thermal photons, otherwise, 
if string-tube model of thermalization takes place, then yield of thermal photons in hard part will be much smaller. 

We estimated the ratio of Direct/Decay photons at LHC energy and find it large enough (20 — 30%) to be mea- 
sured by ALICE PHOS setup, whose expected sensitivity to the direct photons is approximately 5% [Q. So that 
thermal photons considerably contribute into total photon yield. If thermal photons from QGP are evaluated without 
accounting of bremsstrahlung and annihilation with scattering on third particle, then yield of thermal photons from 
hadronic gas exceeds yield of thermal photons from QGP l^]. Otherwise, if one takes into account these higher order 
as corrections, then thermal photons from QGP dominates over photons from hadronic gas pq]. 

The authors would like to thank V.I. Manko for interest and useful discussions. This work was supported by grant 
RFFI 96-15-96548. 
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